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Abstract 

We investigate the dynamical information exchange between a two-state system and its environ- 
ment which is measured by von Neumann entropy. It is found that in the underdamping regime, 
the entropy dynamics exhibits an extremely non-Markovian oscillation-hump feature, in which 
oscillations manifest quantum coherence and a hump of envelop demonstrates temporal memory 
of bath. It indicates that the process of entropy exchange is bidirectional. When the coupling 
strength increases a certain threshold, the hump along with ripple disappears, which is indicative 
of the coherent-incoherent dynamical crossover. The long-time limit of entropy evolution reaches 
the ground state value which agrees with that of numerical renormalization group. 
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The non-equilibrium evolution of open quantum systems is one of the most challenging 
and intriguing problems of contemporary research in both theoretical and experimental 
physics. The transient dynamics can be harnessed and controlled desirably in quantum 
information processing. The first step to manipulate it is to understand how it evolves in a 
short time interval. It is known that the correlations of open system with its surrounding 
environment lead to finite lifetime of quantum superpositions, which give rise to the evolution 
from pure states into mixed ones. It is often stated that decoherence causes the system 
to become entangled with its environment, and the entanglement between them can be 
measured quantitatively by von Neumann entropy l| . The main questions which now arise 
are: How does the entropy or quantum information flow from system to environment? 
After the state of system is initialized as a pure state without entanglement, how does the 
entropy evolve to its long-time limit (monotonously or not)? Is the process of information 
transfer between bath and system, unidirectional or bidirectional? In this paper, as far as 
we know, it is the first time to show the time evolution of entropy for open system which 
exhibits extremely non-Markovian characters and point out the process of entropy exchange 
is bidirectional in the underdamping region. 

The dissipative two-state system (TSS), which is also called the spin-boson model, as a 
simple paradigm of open system, is a generic model which can be widely used to describe a 
large number of physical and chemical processes, such as the defect-tunneling and electron 
transfer, and applied to clarify very interesting quantum phenomena, such as decoherence 
and dephasingj2|, |3|. The open system inevitably encounters decoherence which renders a 
quantum superposition state to decay into a classical, statistical mixture of states. Thus, 
derivation of the reduced density matrix is a central goal in order to describe its evolu- 
tion. Based on weak coupling assumption, a Markovian master equation could give its 
dynamics. However, strong interactions with low-temperature reservoirs give rise to large 
system-environment correlations which generally result in a failure of the Markovian ap- 
proximation. In this case, the system dynamics possesses long memory times and exhibits 
a pronounced non-Markovian behaviors j4|. Thus, it is significant to show the temporary 
evolution of system by a non-Markovian approach, especially, in the case of the strong cou- 
pling to its bath. The non-Markovian approach can investigate not only more complicated 
situation where Markovian approximation is unreachable but also different spectral densities 
between the system and the environment. 
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The entanglement entropy has been considered in many works. They mainly focus on the 
static or ground state properties and find some important results js-?]. Costi and McKenzie 
used a numerical renormalization group (NRG) treatment to study the entropy of the ground 
state as a function of couplingisi]. Recently, Hur and coworkers applied the NRG to study 
the quantum phase transition and found that there is a cusp in the entanglement entropy 
accompanying with quantum phase transition jol. As explicitly pointed out by Costi and 
McKenzie, it is fascinating to show how entropy varies with time after the qubit is initially 
prepared in a certain state without entanglement |5|. In order to calculate the dynamics 
of entanglement entropy we present an analytical approach based on a unitary transforma- 
tion method without the Markovian approximation. It works well in the parameter regime 
< a < 1 and < A < Uc, and can reproduce well known non-perturbation results obtained 
by various methods, such as the coherence- incoherence transition, which have been studied 
in our previous work [8]. The approach does not invoke the rotating- wave approximation so 
as to take into account the effects of count-rotating terms on transient dynamics. Since the 
quantum manipulation can be effectively made in the coherent region, we would give the 
evolution of entropy for this regime at T = 0. We find the non-Markovian entropy evolution 
with pronounced small oscillations feature in the weak coupling, which demonstrates quan- 
tum coherence. As the coupling increases, a hump along with ripple clearly emerges in the 
short time characterizing the temporal memory of bath, eventually the dissipative effects 
quench some oscillations and the hump near the coherence-incoherence transition. 

The spin-boson model reads[2, |3| 

H = -yCT^ + Ukblbk + \Y1 3kibl + bk)(Tz- (1) 

k k 

Standard notations are used Sj], A is the bare tunneling matrix and Qk the coupling constant. 
The effect of environment is determined by its spectral density: gl6{uj — ijJk) = 2aijj9{uJc — 
co), where a is the dimensionless coupling constant, Uc is a cutoff frequency and 9{x) is 
the usual step function (In the work the spectrum is Ohmic type, and we set h = = 
1). Although the model seems quite simple, it is in general not exactly solvable and a 
large variety of approximate analytical and numerical methods have been proposed and 
implemented to study its ground state and dynamics jsl- jisl] . 

Unitary transformation. A unitary transformation, which is defined as H' = 
exp{S)H exp{—S), is applied to H in order to take into account the correlation between 



the spin and bosons js, l9|. The form of generator is proposed, 

Qk 



S = ^^jM-bk)^., (2) 



where a fc-dependent function is introduced 
the end and the result is H' = Hq + H[ + H2, 



8". The transformation can be performed to 



^0 = + Y.^kblb, - - 6), (3) 

k k 

H[ = ^ - + ^'^)^^ - ^'^yB, (4) 

k 

^2 = - (cosh{5} -T])- ^lOy (sinh{fi} - ^B) , (5) 



where B = ^^kibl — bk) and A^. = 77 A. The renormahzied factor of tunnehng is 

2 

1] = exp[—J2kT^^k\- ^0 is the unperturbed part of H' and, obviously, it can be solved 

k 

exactly since the spin and bosons are decoupled. The ground state of Hq is \go) = \si)\{Ok}) 
(o'xlsi) = |si),|{Ofc}) is the vacuum state for every boson mode Uk = 0). H[ and H2 are 
treated as perturbation and they should be as small as possible. For this purpose rj is 
determined to make Ttb{pbH2) = 0, where pb is the density operator of bath. Besides, 
is determined as 

UJk + 

and because of this form H[ is rewritten as 

H[ = J2Vk[blcr^ + bka^], (7) 

k 

where Vk = ^rQkik/^k and cr_ = (ctz — iay)/2, (T+ = (cr^ + When T = it is easy to 

check that B'^goj = 0. This is essential in our approach. 

In our treatment Hq is treated as the unperturbed Hamiltonian, in which the tunneling 
has been already renormalized by r] coming from the contribution of diagonal transition 
of bosons. H[ is the perturbation relating to the non-diagonal transition of single-boson, 
and H2, containing all other multi-boson non-diagonal transitions, is omitted because its 
contribution to physical quantities is 0{gl) and higher. Note that < < 1- C,k measures 
the adiabatic intensity of the particle interacting with its environment [sl . ~ 1 if ^ A^, 
while i^fc ^ 1 for 0;^ ^ A^. In addition, by the choice of 1^^, H[ has taken into account the 



effects of counter-rotating terms. In other words, the bare couphng (?fc/2 in the original 
Hamiltonian is replaced by the renormalized coupling Vk after the unitary transformation. 

Density operator. In order to show the quantum dynamics, we would first give the 
density operator in Schrodinger representation, pssit) with Hamiltonian H, where the sub- 
script SB stands for the spin-boson model. For the transformed Hamiltonian H' the density 
operator is p'ssit) = psB(t)e~^ . The density operator in the interaction representation is 
Psb{^) — 6xp(i-ffo^)p55(t) exp(— if^Qt). By the equation of motion for p'ssi^) [l6j, we obtain 
the master equation 



d 
dt 



pUt) = - f^MH[{t). [H[{t'),p'^sit')PB]]dt'. (8) 
Jo 



where p'/(t) = TiBp'sB^t) and H[{t) = exp{iH'f^t)H[ exp(— zi/gt). It is known that one can 
arrive at the Born-Markov approximation equation neglecting retardation in the integration, 
i.e., Ps{t') is replaced by p'g{t). Our treatment is beyond this approximation. 

/io\ 

At t = 0, the usual initial density operator is Psb{0) = ps{0)pB = Pb- Then 

{oo) 

we can get the initial condition for our calculations: p'sb{^) = p'sb^^) = ^^PsB{0)e ^ leads 
to Ps{0) = ( I . The calculation is up to the second order gl and and the details are 

oy 

/ P'll Pl2 

shown in the Appendix. The solution of reduced density operator p'g{t) = | | is 

P21 P22 

P'nit) - P22W = ^ y e-^-'duF*iuj) + — j e^-'duF{u;), (9) 



P'i2(i) + P2i(i) = 1 - TT- / -X T' (10) 



27ri 



where Fiyi) = (u; — A^. — X]fc u;-^fc-io+ )~^' "^^^ real and imaginary parts of^i^V^/{uj — iO^ — 
Uk) are denoted as 

LO + lAr t -|- Ar iO + A.^ [A^^Wc — UJ) j ) 

j{u) = 2antu (0 < u; < u,), (12) 

respectively. 
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Dynamical quantities. In what follows we calculate the dynamical quantities, 
{(Trit)) = Tt sTt b[ps bo't]{t = X, y, z). The reduced density operator of the original Hamilto- 
nian H is psif) = T^T^BPsBit), which can be expressed as ps(t) = |[1 + '^^{o'r(t))ar]- First, 
we calculate (^^(t)) which is usually denoted as P{t) in the literature. 

Pit) = Tr.TrsM)..) = f . /'"^"t/i y y (13) 

since Ti bPb = 1- The integration in Eq.Q can be done approximately by the residue 
theorem, P(t) = cos(ci;ot) exp(— 7^), where Uq is the solution of equation wq — A,. — i?(a;o) = 0, 
and 7 is the Wigner-Weisskopf approximation of 7(w): 7 = jaA^. The solution is real 
only when a < ac, oic =_(1 + ^r/^c)/'^- It becomes the well-known result Oc = 1/2 in the 
scaling limit A/wc <^ 10, Q]- For a > there is no real solution uq and it means that 
a = ttc determines the critical point corresponding to the coherent-incoherent transition. 
The coherent regime a < ac can be divided into the underdamping part and the overdamping 
one by a criterion ujq > 7(c<;o) (underdamping), or Uq < 7(^0) (overdamping). In the scaling 
limit A/cOc ^ 1 the point where uq = 7(^0) is at a* = 0.325, which is very close to previous 



results a = 1/3 or 0.3 12|, [iSj. From Eq. (1) one can get a relation between {(Jylt)) and 
i'^zi't)), {cyii)) = ~\'it{^z{'t)) 1 since i[H^az] = —Aay. {crx{t)) can be calculated in the 
following, 

where ^(0;) = R{Ar + u) — R{Ar — u) and T^u) = 7(Ar + u) + 7(Ar — u). One can check 
that the initial conditions (crx.(0)) = 0, (cry(O)) = 0, (cr^(O)) = 1 are exactly satisfied. Besides, 
((Ja;(oo)) = T], ((Ty(oo)) = 0, (cr^(oo)) = 0, which are the correct results for thermodynamic 
equilibrium state [5|. 

Entropy of entanglement. The entropy (indeterminacy of the state) is a measure 
of the missing information compared with the pure state of the composite system. The 
more lost information about the composite state, the more information is contained in the 
correlation between the substates. The greater the entropy of system, the more strongly is 
the pure state of composite system correlated and thus entangledj^]. To see what happens 
to the coherence properties due to the interaction between the system and its surrounding 
starting from a pure state, we use the von Neumann entropy. It is defined as S{t) = 
—Tt{Ps log2 ps), which is a measure of the entanglement between them. It may be expressed 



in terms of the eigenvalues A±(t) = 1/2 ± ^J {(7x{t)y + {(yy{t)y + ((j^(t))2/2 of the density 
operator ps as, S{t) = — A+ logg A+ — A_ logg A_. 

Form the Hamihonian (Eq. 1), it predicts that ((jj,(oo)) = 0, ((^^(oo)) = 0, and only 
{ax{oo)) is nonvanishing in the delocalized phase, which verify our obtained results. So, the 
entropy in long-time limit Seq is given by A± = 1/2 ± (cr^) /2, which is shown in Fig. 1 along 
with the NRG results j^. As a increases, Sgg becomes large. When a — )■ 1, Seq tends to 
one. In the scaling limit, for a > 1, (ax) = and the system remains its initial state, thus 
\{(Jz)\ = 1 and Seq = 0. In other words, the transition between localized and delocalized 
phase occurs at a = 1 and the entropy decreases from unity to zero abruptly. 

In order to calculate the entropy, (cra;(oo)) can also be evaluated using the NRG applied 
to the equivalent anisotropic Kondo model. The NRG data shown in Fig.l are taken from 
Ref. [sl. It is seen that for small tunneling our result is in good agreement with those of NRG. 
However, with increasing large tunneling some discrepancies appear for moderate values of 
the coupling. We think that it comes possibly from the NRG discretization 13j . 

The dynamics of entanglement entropy displays extremely non-Markovian features. Fig- 
ure 2 shows S{t) for different couplings with A = O.luc- For a < a*, the entropy increases 
non-monotonically from zero to a finite value (Seq) with explicit oscillations, and would not 
come close to saturation in the short-time interval which means that quantum coherence is 
not directly destroyed by the bath. At the same time, the envelope of entropy exhibits a 
hump characterizing short-time memory of bath. On the other hand, the oscillation-hump 
feature demonstrates that the process of entropy exchange is bidirectional. To better under- 
stand the nature of the oscillation-hump feature and the large contribution from quantum 
fluctuations, we should consider the elements of the reduced density matrix. {o-y{t)) and 
(czit)) exhibit oscillations which represent coherence. So, the dominant contribution to the 
oscillatory signal comes from {(7y{t))'^ + {a zlt))"^ , while the trend of entropy evolution ascribes 
to {ax{t)). Thus, the oscillation of entropy shows the coherent evolution in coherent regime. 

As coupling increases, oscillations become obviously weaker with small amplitudes and 
the envelop of entropy rises rapidly with small hump (To see Fig. 2b) due to the effects of 
strong dissipation. Near the crossover from coherent to incoherent regime a ~ 1/2, entropy 
shape displays faster rising behaviors without oscillation and the hump disappears, which is 
an important character corresponding to the coherent-incoherent crossover. Note that Seq 
is analytic and continuous at ac because no phase transition happens at this point while 



7 



the dynamical crossover from damped oscillatory to pure decaying behaviors takes place. 
Thus, we can not extract a distinguishable feature of this crossover from Seq because of 
its character of thermodynamical equilibrium even if a > Seq is near to its saturation. 
Therefore, only transitory dynamics of entropy could give the indicator of the crossover even 
though Seq might also be regarded as an interesting order parameter to mark quantum phase 
transitions. 

The evolution of entanglement entropy is very different from that of Markovian approxi- 
mation. The dynamics of S{t) is shown in Fig. 3a with several tunnelings for a = 0.2 as well 
as the corresponding Markovian results. In the Markovian evolution, the system undergoes 
a smooth and fast relaxation to its final statistical mixture. It is found that there is no 
short-time oscillations in the Markovian evolution. Thus, the transient oscillatory behaviors 
of entropy dynamics can not be correctly described by the Markovian approximation. Nev- 
ertheless, in the long-time limit, Markovian results are consistent with Seq as expected. The 
oscillation of the entropy is a hallmark of non-Markovian dynamics in the coherent regime 
which is unexpected in the Markovian dynamics. From the scaled entropy S{t)/Seq in Fig. 
3b, one can see that the entropy displays almost synchronously with different amplitudes of 
oscillations for any tunneling and eventually goes wiggly down to Seq- It indicates that the 
system exchanges entropy frequently with its environment in the short time. In this case, 
the oscillations are more pronounced for the enhancement of coherence involved transition 
between two states. From another point of view, the ability of exchanging information be- 
comes strong for the system with increasing tunneling and it remains coherence for a longer 
time. (Note that the unit of time is rjA, which becomes explicitly larger with increasing 
tunneling.) 

In the coherent regime, a sufficient number of quantum manipulations can be performed 
within the coherent time. The need to maintain quantum coherence during the operation is 
especially difficult to achieve in solid state systems such as quantum dots which couple rela- 
tively strongly to uncontrollable environmental degrees of freedom, leading to decoherence. 
Only in the underdamping regime, the quantum control has more efficiency. The promising 
experimental proposal that entanglement entropy can be measured in Cooper pair box or 
quantum dot scheme is suggested by Kopp and Hur recently [Q]. We really expect that ex- 
perimental setup is capable of testing our predictions and such measurements would provide 
a proof of the existence of oscillations in the entropy evolution although it is not easy to 
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probe small signals in the background of noises and thermal fluctuations. 

Summary. The entanglement entropy dynamics of dissipative TSS is studied by means 
of the analytical approach on the basis of a unitary transformation. Analytical results of 
the quantum dynamics, described by the Ps{t), is obtained for the general finite A/uc case. 
The entanglement entropy evolution from a pure state is shown with explicit non-Markovian 
features. Our approach is quite simple and tractable without spectral structure dependence, 
and it could trigger many future applications in other more complicated coupling systems 
with realistic spectrum function, such as superconducting qubit with Lorentz spectrum. 

Here are a few words about the key ingredient of the approach. The purpose of our 
unitary transformation is to find a better way to divide the transformed Hamiltonian into 
unperturbed part Hq, which can be treated exactly, and perturbation ones H[ + H2, which 
may be treated by perturbation theory. In Hq the tunnelling has been already renormalized 
by 77 which comes from the contribution of diagonal transition of bosons. H[ is related to the 
non-diagonal transition of single-boson and all other multi-boson non-diagonal transitions 
are contained in H2. If one treats the coupling term in the original Hamiltonian H as the 
perturbation, the dimensionless expanding parameter is gl/ool- For Ohmic bath s = 1 it is 
2a /u which is logarithmic divergent in the infrared limit. By choosing the form of rj and 
introducing the function in the unitary transformation it is possible to treat H[ and H2 as 
perturbation because of the following reason. On account of the form of 77 H2 can be treated 
as perturbation because its contribution is zero at second order of gk- The effect of the 
coupling term in H' {H[) can be safely treated by perturbation theory because the infrared 
divergence in the original perturbation treatment for H is eliminated by making choice of 
the function form ^fc. The expanding parameter (s = 1) is gf^f/uf ~ 2au/{u + r/A)^, 
which is finite in the infrared limit. Besides, our approach is well checked not only by the 
initial values of the correlation functions and entanglement entropy, such as P{t — 0) — 0, 
S{t = 0) = 0, and their long time limits such as P(oo) = 0, S{oo) — Seq. 

This work was supported by the China National Natural Science Foundation (Grants 
Nos. 10734020 and 90503007). 
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Appendix 



In this Appendix we list the details of solving the master equation (8). The integration 
in Eq.(8) can be done as follows, 



- / Trs[H[{t),[H[it'),pUt')pB]]dt' 
Jo 

-nk(J+Ps{t')a- + {uk + l)ps{t')(J+aJ\ exp[i(wfc - Ar)(t - t')] 
+ [{uk + l)(7+(7_p'i(i') - nfc(7+p'i(i')(7- - [uk + 

+nkp'^{t')a_a+\ exp[-i(a;fe - A,)(i - t')]] , (Al) 

where Uk — l/[exp(/3a;fe) — 1] is the Bose function. Thus, Eq.(8) can be solved by the Laplace 
transformation. If we denote 

( P'li P'li 



PsiP) 



P'ii p'li 



the solution of Eq.(Al) is 



// // _ 1/2 1/2 

,1 _ ^ V2 V2 M3) 



/ ^ cothK/2T) \ 



Using the relation between Schroedinger and interaction representation and making the 

Laplace inverse-transformation, we can get 

Piiit) - P22it) = cos(A,t)(p'A(t) - p'l,{t)) -ism{Art){pUt) - p'L{t)) 



(A5) 
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The integration path is on a hne parallel to the imaginary axis of complex p plane from 
p — — ioo to p = 0"*" + ioo. 
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